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Interplay between orbital and spin degrees of freedom is theoretically studied for the phase transi- 
tion to the spin-singlet state with lattice dimerization in pyroxene titanium oxides y4TiSi2 06 (y4=Na, 
Li). For the quasi one-dimensional spin- 1/2 systems, we derive an effective spin-orbital-lattice cou- 
pled model in the strong correlation limit with explicitly taking account of the t2g orbital degeneracy, 
and investigate the model by numerical simulation as well as the mean-field analysis. We find a 
nontrivial feedback effect between orbital and spin degrees of freedom; as temperature decreases, 
development of antiferromagnetic spin correlations changes the sign of orbital correlations from 
antiferro to ferro type, and finally the ferro-type orbital correlations induce the dimerization and 
the spin-singlet formation. As a result of this interplay, the system undergoes a finite-temperature 
transition to the spin-dimer and orbital-ferro ordered phase concomitant with the Jahn- Teller lattice 
distortion. The numerical results for the magnetic susceptibility show a deviation from the Curie- 
Weiss behavior, and well reproduce the experimental data. The results reveal that the Jahn- Teller 
energy scale is considerably small and the orbital and spin exchange interactions play a decisive role 
in the pyroxene titanium oxides. 

PACS numbers: 75.30.Et, 75.10.Jm, 75.40.Cx 



I. INTRODUCTION 

Orbital degree of freedom has attracted much attention 
since it plays key roles in electronic properties of transi- 
tion metal compoundsii*^ The orbital degree of freedom 
couples with the Jahn- Teller (JT) lattice distortion, and 
in many compounds, the energy scale of the JT stabi- 
lization energy is larger than that of the spin and orbital 
exchange interactions^ A typical example is the mother 
compound of colossal magnetoresistive (CMR) mangan- 
ites, LaMnOs! The JT energy scale is ~ 0.1-leV and the 
spin exchange energy scale is '-^ lOmeV. Consequently, 
the orbital-lattice transition occurs at a much higher 
temperature (~ 800K) than the antiferromagnetic (AF) 
transition temperature (~ 140K)i4 In these systems, the 
orbital-lattice physics is dominant, and the orbital and 
lattice orderings modify effective spin exchange interac- 
tions to lead a magnetic ordering in a secondary effect. 
Moreover, the JT distortion suppresses fluctuation effects 
in orbital and spin degrees of freedom. Hence, a large JT 
coupling masks bare interplay between spin and orbital 
in many real materials. 

Quantum and thermal fluctuations in the competition 
between the spin and orbital exchange interactions may 
yield novel phenomena, and therefore, it is highly desired 
to explore systems in which the genuine spin-orbital in- 
terplay appears explicitly. One of promising candidates 
is the t2g electron system such as titanium and vanadium 
oxides. In the t2g systems, the JT coupling becomes 
smaller than in the eg systems such as CMR mangan- 
ites because spatial shapes of t2g orbitals avoid directions 
of surrounding ligands in the octahedral positions. For 
instance, in vanadium perovskite oxides AYO^ {A=La,, 
Ce), the orbital-lattice transition temperature becomes 
even lower than the AF one;^ and novel interplay be- 
tween spin and orbital is proposed theoretically. §ii 



FIG. 1: (a) Lattice structure of pyroxene oxides ylTiSi206 
(j4=Na, Li) . Chains of TiOg octahedra are separated by Si04 
tetrahedra. (b) Skew edge-sharing chain structure of TiOg 
octahedra. Balls in the center of each octahedron denote Ti 
cations, and oxygen ions are on the corners of the octahedra. 
The octahedra share their edges in the xy and yz planes al- 



orbitals of t2g electrons, and t^^ {t^' 
integrals between d^y (dyz) orbitals. 



) is the a-bond transfer 
See the text for details. 



Pyroxene titanium oxides ATiSi2 06 (A=Na, Li) are 
typical examples of the t2g electron systems where such 
interplay between spin and orbital is expected^ The lat- 
tice structure of these compounds consists of character- 
istic one-dimensional (ID) chains of skew edge-sharing 
TiOg octahedra as shown in Fig.n(a). The TiOg chains 
are bridged and well separated by Si04 tetrahedra, and 
therefore, interchain couplings are considered to be much 
weaker than intrachain interactions. In each TiOg chain, 
as shown in Fig. ^ (b) , the edges of octahedra in the xy 
and yz planes are alternatively shared between neighbor- 
ing octahedra, which leads to the zig-zag structure. Each 
magnetic Ti^+ cation has one d electron in these insulat- 
ing materials. Hence, we may consider that a quasi ID 
spin- 1/2 system is realized. 

Pyroxene titanium oxides show a peculiar phase tran- 
sition. The magnetic susceptibility shows a sharp drop 
at Tc = 210K and 230K in NaTiSi206 and LiTiSi206, 
respectively, which indicates the transition to a spin- 
singlet state with a finite spin gap.^ Below Tc, a dimer- 
ization of Ti-Ti distances along the chain was observed 
by the X-ray scatteringiS These remind us of a spin- 
Peierls transitioniiS However, above Tc, the magnetic 
susceptibility shows an unusual temperature dependence 
which clearly deviates from that of other spin-Peierls 
compounds. In spin-Peierls systems, the dimerization 
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is caused by the magnetoelastic coupling, and therefore, 
the transition occurs when short-range spin correlations 
develop enough to drive the lattice dimerization. This 
development of spin correlations is manifested in a broad 
peak of the magnetic susceptibility above Tc, and the 
peak temperature gives a rough estimate of the spin ex- 
change energy scale. Contrary to this conventional be- 
havior, the magnetic susceptibility of the pyroxene com- 
pounds at high temperatures increases as temperature 
decreases, and suddenly drops at Tc without a clear for- 
mation of the broad peaki^ This suggests a breakdown 
of the simple spin-Peierls picture in these pyroxene com- 
pounds. 

For the peculiar transition to the spin-singlet state, an 
importance of the t2g orbital degree of freedom has been 
pointed outi^ Theoretically, a scenario of the orbital- 
driven spin-singlet formation has been explored in spin- 
orbital coupled models for t2g electron systems»ii It was 
proposed that the orbital ordering may modify effective 
spin exchange interactions and induce the spin-singlet 
formation. However, models were considered only for 
systems with corner-sharing octahedra, and hence, it is 
unclear that the argument is applicable to the present py- 
roxene systems with the edge-sharing octahedra. A simi- 
lar scenario has been proposed for the present systems )i2, 
however, the analysis was heuristic and not sufficient to 
conclude the mechanism of the finite-temperature tran- 
sition. In order to clarify the nature of the transition 
and the low-temperature phase, we need more elaborate 
analysis. In particular, it is indispensable to investigate 
thermodynamic properties on the basis of a microscopic 
model. 

In the present study, we will theoretically explore the 
mechanism of the unusual phase transition in the pyrox- 
ene titanium oxides ATiSi2 06 explicitly taking account 
of the t2g orbital degree of freedom. We derive an effec- 
tive spin-orbital-lattice coupled model in the strong cor- 
relation limit, and investigate thermodynamics as well 
as the ground-state properties of the model applying the 
numerical quantum transfer matrix method and mean- 
field-type arguments. As a result, we find that interesting 
interplay between orbital and spin degrees of freedom oc- 
curs in the system, and the interplay gives a comprehen- 
sive understanding of the peculiar properties of the py- 
roxene compounds: Although both orbital and spin cor- 
relations are antiferro type and compete with each other 
at high temperatures, development of AF spin correla- 
tions with decreasing temperature yields a sign change of 
orbital correlations from antiferro to ferro type. After the 
sign change, the ferro-type orbital correlations grow with 
the antiferro-type spin correlations, and finally, induce 
the spin-singlet formation with a dimerization. Further- 
more, the nontrivial temperature dependence of orbital 
correlations modifies an effective magnetic coupling, and 
results in a non Curie- Weiss behavior of the magnetic 
susceptibility. We show that our model with realistic 
parameters reproduces the peculiar temperature depen- 
dence of the magnetic susceptibility in experiment. 



The organization of this paper is as follows. In the fol- 
lowing section^ applying the strong-coupling approach, 
we derive an effective spin-orbital-lattice coupled model 
for the t2g pyroxene compounds ATiSi2 06. In Sec. IIIII 
we discuss properties of the system in the ground state 
and in the high-temperature limit using mean-field type 
arguments. In Sec. II VI we study thermodynamic proper- 
ties of the effective model by numerical simulations. We 
describe the method in Sec. IIV Al Section IIVBI shows 
the numerical results including quantitative comparisons 
with the experimental data. Finally, Sectionals devoted 
to summary and concluding remarks. 



II. MODEL HAMILTONIAN 

In this section, we derive an effective spin-orbital- 
lattice coupled model for the pyroxene oxides ATiSi2 06, 
whose Hamiltonian reads 



TL = Hs 



JT 



■Ha 



(1) 



The first term describes the intersite exchange interac- 
tions in spin and orbital degrees of freedom, and the 
second term includes the Jahn- Teller type orbital-lattice 
coupling. These two are defined within each ID chain. 
The third term describes the interchain coupling. 

The spin-orbital exchange Hamiltonian Tiso is derived 
from a t2g multiorbital Hubbard model by using the per- 
turbation in the strong correlation limit The skew 
structure shown in Fig. (b) distorts TiOg octahedra in 
the form that four Ti-0 bonds in which the oxygen ions 
are shared with neighboring TiOg octahedra are longer 
than the rest two Ti-0 bonds in each octahedron.^. The 
distortion leads to the splitting of threefold t2g levels into 
a low-lying doublet with dxy and dyz orbitals and a sin- 
gle dzx level when we take the coordinates as shown in 
Fig. n (b). The splitting is estimated as ~ 300meV in 
a pyroxene vanadium oxide which has a similar lattice 
structure.^'* Because of this level splitting, it is reason- 
able to start from the ID Hubbard model with twofold 
degeneracy of d^y and dyz orbitals with tracing out the 
higher dzx level. The Hamiltonian is given in the form 

^Hub = X]mi(^S'^^J"rCj73r +H.C.) 

i,j a, [3 T 



i af3^a'f3' tt' 



(2) 



where i,j are site indices within the ID chain, r, t' are 
spin indices, and a, P — 1 (dxy) and 2 (dyz) are orbital in- 
dices. The first term describes the electron hopping, and 
the second term denotes the onsite Coulomb interactions 
where we use the standard parametrizations, 

Uaf3.a't3' = U'Saa'SpfJ' + Jh ('5a/3"5/3a' + SafsSa'P'), (3) 

U = Un,n = U22,22 = U' + 2Ju. (4) 
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Here, we neglect the relativistic spin-orbit coupling. 

The perturbation calculation is performed in the 
strong correlation limit t"-^ ^ U,U', Jh by taking an 
atomic eigenstate with one electron at each site. In the 
edge-sharing configuration as shown in Fig. ^ (b), the 
most relevant contribution in the transfer integrals t"j^ 
is the overlap between the nearest-neighbor (NN) pairs 
with the same orbitals lying in the same plane, which 
is called the a bond. The tr-bond transfer integrals are 

i+i for NN pairs in the xy plane and tf1_^^^ for those in 
the yz plane as shown in Fig. (b). These two types of 
the transfer integrals take the same value, and we denote 
them by in the following. Other transfer integrals are 
much smaller than to-; in particular, i^^+i — tf\^i = 
between NN sites from the symmetry. In the present 
study, we take account of the u-bond contributions only, 
and neglect other transfer integrals The approxima- 
tion is known to give reasonable results for spinel oxides 
which also have edge-sharing network of octahedraii^ 

The second order perturbation in t„ gives the effective 
spin-orbital Hamiltonian in the form 



7 oAF _ 



(5) 



. oF 
"-1,1+1 



{A + BS,-S,+i)i^-^2T,T,+,j, (6) 

(7) 



where Si is the S = 1/2 spin operator at site i and Ti is 
the Ising isospin which describes the orbital state at site 
z as Ti = -|-l/2(— 1/2) when the dxy {dyz) orbital is occu- 
pied. Note that the cr-bond transfer integral t„, which is 
orbital diagonal and does not mix different orbitals, leads 
to the Ising nature of the orbital isospin interaction. The 
signs in Eq. iQ take + (— ) for the bonds within the xy 
{yz) plane. The coupling constants in Eqs. ©-0 are 
given by parameters in Eq. Q as 
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(8) 
(9) 

(10) 

(11) 
(12) 



where we use Eq. Q. The realistic value of 77 will be 
estimated as ry ~ 0.1 later. Therefore, we consider that 
A, i3, and C are all positive in the following. 

We note that a part of the spin-orbital interactions, 
i.e., in Eq. iQ takes a similar form to the model 



proposed in Ref. ^3 Our model derived from the mul- 
tiorbital Hubbard model contains both ferro and anti- 
fcrro types of orbital interactions as well as the nontriv- 
ial contributions which are missed in the previous model 
in Ref. 12. We will show that these factors play impor- 
tant roles in the thermodynamics. We also note that 
a spin-orbital model similar to Eq. was studied in 
Ref. IllL The model was derived for the corner-sharing 
configuration of the octahedra, while our model is for 
the edge-sharing configuration. 

The orbital-lattice coupling term TYjt in Eq. is 
given in the form 



(13) 



The first term describes the JT coupling where 7 is the 
electron-lattice coupling constant and Qi is the ampli- 
tude of the JT distortion which couples to the remaining 
two- fold degeneracy of the d^y and dyz orbitals. The sec- 
ond term denotes the elastic energy of the JT distortion. 
For simplicity, we neglect the kinetic energy of phonons 
and regard Qi as a classical variable. Here, we note that 
besides the onsite term there may be intersite interac- 
tions of JT distortions such as 'YliijQiQj- However, in 
the self-consistent scheme described in Sec. IIV Al which 
we will employ in the present numerical study, the inter- 
site interactions do not affect the results except for a shift 
of the critical temperature. Such effect can be renormal- 
ized into the 7 term in Eq. H13|) . and therefore, we do 
not explicitly include the intersite term in the present 
Hamiltonian. 

In addition to the above two terms Tiso and Tijx, we 
also consider the interchain coupling term T-L± in Eq. 
We may consider two contributions to H.^. One is the 
spin-orbital exchange interaction arising from the inter- 
chain transfer integrals of electrons, and the other is the 
interchain interaction of JT distortions. The former is 
expected to be negligibly small due to a rapid decay and 
a large spatial anisotropy of the t2g wave functions. We 
therefore ignore it and take account of the latter JT con- 
tribution only. The explicit form of the JT coupling de- 
pends on the details of the lattice structure among ID 
chains, and may be very complicated. In the present 
study, we assume the following simple form 



(14) 



where the summation is taken over the NN sites in four 
neighboring chains. Since the TiOg chains are well sep- 
arated by Si04 tctrahedra in the pyroxene compounds, 
the coupling constant A is expected to be small. In the 
following numerical calculations, we will treat the inter- 
chain JT coupling in a mean-field approximation, which 
is well justified for weakly-coupled ID systems.^* 

Finally, we discuss the realistic values of the param- 
eters. The (T-bond transfer integral is estimated as 
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io- — 0.3eV for pyroxene vanadium oxidesii* Typi- 
cal estimates for Coulomb interaction parameters are 
[/ - 5-6eV and Jh - 0.6 -0.7eVili Hence, rj in Eq. ^ 
is a small parameter of the order of 0.1. From the esti- 
mates of to- and U, the exchange energy scale J in Eq. (jSJ 
is estimated as J ~ 200K. As for the JT couplings 7 and 
A, since there is no experimental estimates to our knowl- 
edge, we treat them as parameters which are determined 
by the comparison between the numerical results and the 
experimental data in Sec. IIV B 21 

III. MEAN-FIELD ANALYSIS 

Before going into the numerical study of the thermody- 
namics of the model Q , we apply mean- field arguments 
to capture the behavior in the ground state and in the 
high-temperature phase. For simplicity, we consider only 
the spin-orbital part Tiso in this section. 

A. Ground state 

In the ground state, we here consider four different 
types of ordered states as shown in Fig. |21 schemati- 
cally; (a) the spin-ferro and orbital- ferro (sF-oF), (b) the 
spin-ferro and orbital- antiferro (sF-oAF), (c) the spin- 
antiferro and orbital-ferro, and (d) the spin-antiferro 
and orbital-antiferro (sAF-oAF) states. For simplic- 
ity, we consider the fully polarized states for the spin- 
ferro, the orbital-ferro, and the orbital-antiferro order- 
ings, where quantum fluctuations do not play a crucial 
role. (Note that the orbital isospin is the Ising spin in 
the present model.) The important point is that as easily 
shown by Eqs. © and Q, the orbital-ferro ordering, i.e., 
(r,r,+i) = 1/4 and (T,) = 1/2 (or -1/2) for all i, discon- 
nects every other bonds. Hence, when the spin coupling 
is antiferromagnetic, the orbital-ferro ordering leads to 
the spin-singlet formation, i.e., {Si ■ Si+i) = —3/4 for re- 
maining isolated NN pairs. Therefore, the ordered state 
(c) is denoted as the spin-dimer and orbital-ferro (sD-oF) 
state. We note that a similar mechanism of the singlet 
formation driven by orbital ordering has been proposed 
for a related modeliJA 

The ground-state energy for each ordered state is cal- 
culated by replacing the spin and orbital-isospin opera- 
tors in Eq. Q by the following expectation values; 

{Si ■ Si+i) = i for sF pairs, (15) 

{St ■ Si+i) = for sD pairs, (16) 

{Si ■ Si+i) = -s for sAF pairs, (17) 

where s is a positive parameter less than 3/4 (we do not 
need the precise value of s), and 

(T,T,+i) = i, {T,) = ±i for oF pairs, (18) 



(a) sF-oF 




(b) sF-oAF 




(c) sD-oF 




(d) sAF-oAF 




FIG. 2: Spin and orbital ordering patterns in (a) the spin- 
ferro and orbital-ferro (sF-oF), (b) the spin-ferro and orbital- 
antiferro (sF-oAF), (c) the spin-dimer and orbital-ferro (sD- 
oF), and (d) the spin-antiferro and orbital-antiferro (sAF- 
oAF) states. Circles denote the lattice sites, and two sep- 
arated spaces inside them denote two different orbital states 
dxy and dy^ as indicated in (a). Arrows denote spins, and 
gray ovals in (c) represent the spin-singlet pairs. Crosses in 
(a) and (c) denote the disconnected bonds by the ferro-type 
orbital ordering. See the text for details. 

(r,T,+i) - (T,) = ±-^^ for oAF pairs,(19) 

respectively. The obtained ground-state energies per site 
for the states (a)-(d) are 

SsF-oF = 0, (20) 

E,F_oAF ^-j(^A+^y (21) 
C 

-£^sD-oF = ~'^'2' ^^^^ 

SsAF-oAF = -J(A-sS), (23) 

respectively. 

We compare Eqs. (|5ni)-(|131l by using Eqs. ©-(HH) and 
obtain 

-oAF or E'sD-of) < -E-sAF-oAF < -E-sF-oF- (24) 
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Hence, the ground state is either (b) sF-oAF or (c) sD- 
oF state. From the equation of i?sF-oAF = £'sD-oF, we 
obtain the critical value of for the transition between 
the two phases as 



1 



(V73 



0.18. 



(25) 



Namely, we have the ground-state phase transition be- 
tween the sD-oF and sF-oAF phases by changing 77; the 
sD-oF (sF-oAF) phase is stable for 77 < rjc [r] > rjc)- 
Hence, for a realistic value rj ^ 0.1, the ground state 
is predicted to be the sD-oF phase within the mean-field 
argument. This will be confirmed by the numerical cal- 
culations in Sec. IIVI We will also show that the sF-oAF 
state is realized for 77 > 77c in Appendix B. 

The 77-controlled phase transition is understood by the 
competition between the spin superexchangc interaction 
and the Hund's-rule coupling. The former comes from 
the perturbation process within the same orbitals and 
favors the spin-singlet stated-. The latter enhances the 
energy gain from the perturbation process between differ- 
ent orbitals and stabilizes the sF-oAF state. It is known 
that the sF-oAF state is favored by a finite Hund's-rule 
coupling (77 > 0) in multiorbital systems with transfer 
integrals = t"^^ ^ between all the NN sitesi^a*2i In 
the present systems, the specific form of the transfer in- 
tegrals due to the zig-zag lattice structure gives a chance 
to stabilize the sD-oF state in the small r] regime. 



B. High-temperature para phase 

At high temperatures, both spin and orbital are disor- 
dered. We here examine spin and orbital correlations in 
the para phase by a mean- field type argument. 

First, to focus on the spin degree of freedom, we replace 
the orbital isospin operators in the effective model Q 
with their mean values, i.e., Ti (Ti) = and TiTi^i — > 
(TiTi+i). The resultant effective spin Hamiltonian reads 



where 



P —J 



K: 



-{C-2B) + {C + 2B){T,T,+i) 



-{C + 8A) + iC-8A){T,T,+,) 



(26) 

,(27) 
.(28) 



In the high-temperature limit, {TiTi-^-i) becomes zero and 
the effective spin exchange constant becomes 



(29) 



which is positive for 77 ~ 0.1 and independent of i. Hence, 
we end up with a ID AF spin Heisenberg model. 



On the other hand, we can also consider a reduced 
Hamiltonian for the orbital degree of freedom by replac- 
ing the spin operators in the model © with their mean 
values. The result is 

(30) 



where 
T° 



K° 



= J 



2A- 



C 



(C + 2B)(5,-5.+i) 



^\[2A+^]-{C-2B){S,-S,+{) 



(31) 
(32) 

(33) 



In the high-temperature limit, {Si ■ Si+i) = and the 
effective isospin coupling constant becomes 



J°,+i(T^oo) = J(2A-5 



(34) 



This coupHng constant is also positive for -q ^ 0.1, and 
therefore, we obtain a ID antiferro-type Ising isospin 
model. [Note that the second term in Eq. (|30|l cancels 
out when {Si ■ Si+i) = 0.] 

The above arguments indicate that both spin and 
orbital correlations are antiferro type in the high- 
temperature limit in our effective model On the 
contrary, as discussed in Sec. lIII Al the ground state is ei- 
ther the sF-oAF or sD-oF state. This suggests that either 
spin or orbital correlation has to change from antiferro 
to ferro type as decreasing temperature. It is implied 
by Eq. H27|) that development of antiferro-type orbital 
correlations {TiTi+i) < may change the sign of Jii^i 
from positive (antiferro) to negative (ferro). In a similar 
way, Eq. suggests that development of AF spin cor- 
relations {Si ■ Si+i) < may induce the sign change of 
Jii^i- In this manner, the spin and orbital correlations 
compete with each other. The mean-field-level argument 
is clearly insufficient to clarify the competition, and we 
will employ the more powerful numerical analysis in the 
next section. 



IV. NUMERICAL ANALYSIS 
A. Method 

To study thermodynamic properties of the model 
we apply the quantum transfer matrix (QTM) method^ 
combined with the mean-field treatment of the JT dis- 
tortions. Here, we describe the scheme of our analysis. 

Since the present model Q is highly ID anisotropic, 
it is justified to treat the weak interchain coupling Ti.^ as 
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the mean fieldi^ in the form 

7i± = A^Q4j=-^|A|^g?, (35) 

where Qj is the mean-field value of the JT distortion at 
site j and z is the number of NN chains, i.e., z = 4 in 
the present materials. As a result, the total Hamiltonian, 
7iso+'^jT+'7^j_, is reduced to an effective ID spin-orbital- 
lattice coupled model in the form 

i i 
i i 

whe re the JT coupling and distortion are rescaled as 7 = 
7/-\/l — 22;|A| and Qi — ^1 — 2z\\\Qi, respectively. For 
later use, we define the JT stabilization energy as 

AjT = 7V8, (37) 

which is the energy gain by the JT distortion in the ab- 
sence of Tiso- 

The optimal values of the JT distortions {Qi] are de- 
termined in a self-consistent way. We start from an initial 
guess of {Qi} and calculate the expectation values of the 
isospin {Ti) by applying the QTM method to the effective 
ID model Ti,. Note that the QTM calculation is numer- 
ically exact and includes all the fluctuations in spin and 
orbital degrees of freedom for a given set of {Qi}. See 
Appendix A for the details of the QTM method. The 
obtained values of (Ti) are used to determine {Qi} self- 
consistently. The self-consistent equation is obtained by 
the energy minimization 5{TL)/5Qi = which gives 

Qr" = -7(T.). (38) 

The values of {Q^i™} are used as inputs for the next 
QTM calculation. We iterate the procedure until {Qi] 
converge to optimal values. 

Thermodynamic properties are obtained for the opti- 
mal values of {Qi{. We calculate the magnetization per 
site, 

m = \Y.^St), (39) 

i 

and the average values of NN spin and orbital correla- 
tions, 

= \ Y.^S,-S,+i), (40) 

i 

Co = ^^(r,r,+i), (41) 

i 

respectively, where L is the length of the chain. We also 
define the NN spin correlations on odd and even bonds 



as 

iGcvcn 

respectively, to detect the spin dimerization. See 
Eqs. (|A6I) - (|A10() in Appendix A for the calculation of 
these quantities. Besides, to calculate the magnetic sus- 
ceptibility X, we perform the calculation within the same 
framework for the system under the external magnetic 
field whose Hamiltonian reads 

H-giiBHY,S!, (44) 

i 

where g is the g-factor and /is is the Bohr magneton, 
and H is the external magnetic field. The susceptibility 
is obtained from the numerical derivative of m as 

^ m(AH)-m(0) 

where we use AH = 0.001J/((7/is). 

The QTM method allows us to treat the system in the 
thermodynamic limit i — > 00 directly (see Appendix A). 
Instead, the results contain a systematic error due to a 
finite Trotter number M. We need to check carefully the 
convergence of the results with increasing M . The values 
of AI used in the following calculations are up to Af = 4 
for the case of 7 = while up to M = 3 for 7 > 0. 
Although these values of M appear to be rather small, 
we will show that the A/-convergence is excellent for the 
temperature range in the present study, T <; O.IJ. 

B. Results 

In this section, we show our numerical results for ther- 
modynamic properties of the model Q obtained by the 
method in the previous section. In Sec. lIVBTl we clarify 
nontrivial feedback effects between orbital and spin de- 
grees of freedom in the absence of the JT coupling. In 
Sec. IIV B 21 we switch on the JT coupling, and compare 
our numerical results with the experimental data. 

1. Interplay between orbital and spin 

First, we show our results in the absence of the JT 
coupling (7 — 0). In this case, the model is reduced to 
the ID spin-orbital model Hso [Eq. ©], and we can ob- 
tain thermodynamic properties by performing the QTM 
calculation without any self-consistent iteration. 

Figure O shows the results of the magnetic susceptibil- 
ity X in Eq. H45|) for several values of 77. Note that the re- 
sults for different values of the Trotter number M = 2, 3, 
and 4 almost coincide with each other as shown in the 
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FIG. 3: Temperature dependence of the magnetic suscepti- 
bility X for r] — 0.0,0.05,0.10, and 0.12. Symbols represent 
the results for M = 4 while solid and dotted curves are those 
for M = 3 and M = 2, respectively. 



figure, which ensures the convergence of the QTM re- 
sults in this parameter range. The magnetic suscepti- 
bility at high temperatures increases as temperature de- 
creases, and exhibits a broad peak at some temperature. 
At lower temperatures, it decays rapidly and approaches 
zero as r ^ 0. While the peak becomes sharper and 
shifts to a lower temperature region as rj increases, the 
rapid decay of x low temperatures is commonly ob- 
served. This indicates that our spin-orbital model Tigo 
in Eq. jSJ exhibits the spin-singlet ground state in the 
small rj regime as predicted in Sec. lIII.^ (See Appendix 
B for the larger rj regime.) 

To clarify the nature of the system in more detail, in 
Fig. 21 we show the results of NN spin and orbital corre- 
lations Cs and Co defined in Eqs. (|4U|I and (|41|l . respec- 
tively. The M-dependence of the results is negligible also 
for these quantities. We note that at finite temperatures 
no true long-range order appears in the ID spin-orbital 
model Hso, and hence the NN spin correlations are uni- 
form, i.e., Cs = Cf'^ = C|™". As T ^ 0, the spin and 
orbital correlations converge to Cs — — 3/8andCo — 1/4, 
respectively. These are indeed the values expected in the 
sD-oF ground state shown in Fig.|21(c), where {Si ■ Si+i) 
takes the value of —3/4 or alternatively from bond to 
bond while (Tf T^^j^) = 1/4 for all bonds. Hence, we con- 
clude that the ground state of the system is the sD-oF 
state in Fig.|21(c) for the realistic values oi ij ^ 0.1, and 
the rapid decay of x at low temperatures in Fig.|21is due 
to the spin-singlet formation in the ground state. 

Let us examine the finite-temperature behavior of Cs 
and Co in Fig.0] At high temperatures, both the spin and 
orbital correlations are negative, i.e., antiferro type, be- 
ing consistent with the mean-field prediction in Sec. lIIIBl 
Here, the spin and orbital antiferro correlations compete 
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FIG. 4: Temperature dependences of (a) the nearest-neighbor 
spin correlation Cs and (b) the nearest-neighbor orbital corre- 
lation Co for ri = 0.0, 0.05, 0.10, and 0.12. Symbols represent 
the results for M = 4 while solid and dotted curves are those 
for AI — 3 and M — 2, respectively. 



with each other since the spin- antiferro correlation fa- 
vors the orbital-ferro one and vice versa as indicated in 
the form of Hso in Eq. ((Sj. We find that in this range of 
rj the spin-antiferro correlation grows more rapidly than 
the orbital correlation as T decreases. This growth of Cs 
suppresses Co and causes the sign change of Co from an- 
tiferro to ferro type. Once Co becomes ferro type, Cg and 
Co develop cooperatively, and eventually the ferro-type 
orbital correlation induces the spin-dimer formation in 
the ground state. These results clearly indicate that our 
model Hso shows a nontrivial feedback effect between or- 
bital and spin degrees of freedom at finite temperatures. 

The strong interplay between orbital and spin also 
shows up in the temperature dependence of the magnetic 
susceptibility at high temperatures. There, we expect a 
deviation from the Curie- Weiss behavior due to the inter- 
play since Jf^^i depends on the NN orbital correlation 
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FIG. 5: Effective spin exchange coupling Jf^^+i in Eq. H27|l 
for r; = 0.0,0.05,0.10, and 0.12. The results are scaled by 
their values in the high-temperature limit, Ji^i+i{T oo). 



(TiTi+i) as mentioned in Sec. 111151 The deviation can 
be monitored by the temperature dependence of the ef- 
fective spin exchange coupling in Eq. (|27|l . since 
it gives an effective Curie- Weiss temperature Qcw- (In 
the mean- field approximation, 8cw — Jfi+i/'^ in the ID 
spin-1/2 system with NN interaction.) 

In Fig. |31 we show the temperature dependence of 
Jfi^i in Eq. I|27|) scaled by the value in the high- 
temperature limit. The results clearly show that 
is temperature-dependent even in the temperature range 
where the Curie- Weiss fittingto the experimental data 
has been attempted in Ref. 13- (The energy scale of J 
will be estimated as 200 — 300K from the fitting in the 
next section.) We need further careful considerations in 
the fitting of the experimental data of the magnetic sus- 
ceptibility to estimate the model parameters. 



2. Comparison with experimental results 

Next, we show our results in the presence of the JT cou- 
pling, i.e., for the case of 7 > in Eq. (|36|) . and compare 
them with experimental data quantitatively. The nu- 
merical results shown in the following are obtained from 
the self-consistent iteration of the QTM calculation with 
M = 3, and we have confirmed the Af -convergence of the 
results. 

Figure IHI shows the results of the magnetic susceptibil- 
ity X. The results show a singularity at some tempera- 
ture and a sudden drop below there, which corresponds to 
the phase transition to the low-temperature sD-oF phase 
as discussed later. For comparison, we plot the exper- 
imental data for NaTiSi2 06,^^ which increase as T de- 



X (emu/mol) 




FIG. 6: Comparison between numerical and ex- 

perimental results of the magnetic susceptibility x- 
Solid, dotted, and dashed curves represent the numer- 
ical results for ( J, 7?, p, Ajt) = (213K, 0.05, 1.85, 86K), 
(250K,0.10, 1.87, 90K), and (296K, 0.12, 1.90, 89K), respec- 
tively. Bold gray curve represents the experimental data in 
Ref. y. 

creases at high temperatures and exhibit a sharp drop at 
Tc — 210K. The numerical results are obtained for several 
typical parameter sets estimated by the following fitting 
procedure: In our model, there are four parameters, i.e., 
the effective exchange coupling J, the ratio of Coulomb 
interactions 77 = Ja/U, the g-factor and the JT cou- 
pling parameter 7. For a certain value of 77, we perform 
the two-parameter fitting by using J and g in the high- 
temperature regime of 400K < T < 700K, where the JT 
coupling 7 is irrelevant in the present scheme of the calcu- 
lation. Using the optimal values of J and g, we determine 
7 to reproduce the critical temperature Tc = 210K. We 
thereby obtain the optimal set of J, g, and 7 for a certain 
value of r/. 

In Fig. El we show the results of the fitting for i] = 
0.05,0.10, and 0.12 as typical examples. The estimates 
of the parameters are (J, Ajx) = (213K, 1.85, 86K), 
(250K,1.87,90K), and (296K, 1.90, 89K) for t] ^ 0.05, 
0.10, and 0.12, respectively. As shown in the figure, the 
numerical results well reproduce the experimental data, 
except for a small deviation near Tc (discussed below). 
Although we cannot determine the best set of the param- 
eter values only from the present fitting, we find that the 
estimates are quite reasonable in this t2g compound: The 
estimates of J ~ 200 — 300K are comparable to the esti- 
mate based on microscopic parameters in Sec. ^ As for 
the g-factor, although there is no estimate for the present 
compound as far as we know, it is known that g becomes 
slightly smaller than two in some t2g compounds)^^*^^ and 
therefore we believe that our estimates of g are reason- 
able. Furthermore, the estimates of Ajx ^ 90K, which 
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are smaller than the exchange energy J, also appear to be 
reasonable for electron systems. (We will comment on 
the magnitude of Ajt in the end of this section.) There- 
fore, we believe that our effective model with realistic 
parameters can describe successfully the peculiar proper- 
ties of the pyroxene compound NaTiSi206. 

We note that the small deviation near Tc can be at- 
tributed to the approximation in the treatment of the JT 
distortion Q^: Our method to determine Qi is not able to 
include effects of the short-range correlations and ther- 
mal fluctuations of lattice distortions, which tend to en- 
hance the spin-singlet fluctuation and suppress the mag- 
netic susceptibility x- It is therefore reasonable that our 
method overestimates x around the critical point, where 
the effects become significant. We believe that better 
agreement near Tc may be obtained by including such ef- 
fects, however this interesting problem is left for further 
study. 

As shown in Fig. El the magnetic susceptibility of our 
model shows an exponential decay at low temperatures 
well below T^. We estimate the spin gap Ag from the 
results below 0.5Tc by the fitting function 



Xcxexp(-A,/r). 



(46) 



The estimates of Ag are 665,640, and 638K for rj = 
0.05,0.10, and 0.12, respectively. In experiment, the spin 
gap As is estimated as ~ 500K^ which is comparable to 
our numerical results. 

Next, we show the results of the spin and orbital cor- 
relations in Fig. [7| Here, we plot only the results for 
the parameter set of (J, 77, Ajt) = (250K, 0.10, 90K) as 
a typical example since essentially the same behavior is 
obtained for other parameter sets in Fig.|HI^^ The corre- 
lations take the same values as those for Ajt = above 
Tc. Below Tc, the NN spin correlations on odd and even 
bonds, C°'^'^ and C"™", take different values and con- 
verge to —3/4 or as T 0. This indicates that the 
translational symmetry is broken in the spin-dimer phase. 
(Note that one of the doubly-degenerate ordered states is 
selected depending on the initial set of {Qi} used in the 
self-consistent calculation.) On the other hand, the NN 
orbital correlation Co is uniform, and goes to the value 
of 1/4 as T — > 0. As shown in the inset of Fig. (b), 
the isospin polarization (T^) becomes finite below Tc, in- 
dicating the long-range orbital ordering. These results 
explicitly show that the system below Tc is in the sD-oF 
phase in Fig. [51(c). 

In Fig. |S1 we plot the effective spin exchange coupling 
Jii+i, corresponding to Eq. lfT7|) in the case of Ajt = 0. 
For a finite Ajt, we have to keep the linear terms of (T;) 
in the reduced spin Hamiltonian in Eq. (|26|) . which give 
rise to an additional term - J(-l)'C/2 ((T^) + (T^+i)) to 
Jli+i in Eq. ^7^. As shown in Fig. (HI the value of Jli+i 
deviates from that for Ajt = below Tc, and takes two 
alternative values from bond to bond. This clearly shows 
the alternating behavior of the magnetic coupling in the 
spin-dimer state. At low temperatures T ^ 150K, the 
dimerization is almost perfect, i.e., Jfi+i on the weak 
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FIG. 7: Temperature dependences of (a) the nearest- 
neighbor spin correlations Q, C°'^'^, and CT™, and (b) 
the nearest-neighbor orbital correlation Co for (J, rj, Ajt) = 
(250K,0.10,90K). Circles show the data for Ajt = for 
comparison. Inset in (b): Polarization of orbital isospin. 



bonds are almost zero. Hence, there the system con- 
sists of almost independent spin-singlet pairs. Note that 
Jfi^i on the strong bonds is surprisingly enhanced up to 
^ lOOOK. The almost perfect dimerization and the en- 
hanced exchange coupling are the unique aspects of the 
present spin-orbital-lattice coupled system. 

Our estimates of the JT stabilization energy Ajt ~ 
90K are considerably smaller than Tc = 210K as well as 
J = 200 — 300K. This suggests that the critical tempera- 
ture Tc is not determined mainly by the JT coupling, and 
that the balance of the JT coupling and the spin-orbital 
intcrsite interaction is important in this t2g electron sys- 
tem. Moreover, we note that the phase transition occurs 
below the temperature where the orbital correlation Co 
changes its sign from negative (antiferro type) to posi- 
tive (ferro type) as shown in Fig. [7\ (b) . This indicates 
an importance of the interplay and the feedback effect be- 
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FIG. 8: Temperature dependence of the effective spin ex- 
change couphng Jf_i+i for (J,r?,AjT) = (250K, 0.10, 90K). 
Circles show the data for Ajt ~ for comparison. 



tween spin and orbital in the present system. Therefore, 
we conclude that the pyroxene compounds are typical t2g 
electron systems where the JT energy scale is relatively 
small and the bare interplay between spin and orbital 
degrees of freedom plays a central role in the thermody- 
namic properties. 



V. SUMMARY AND CONCLUDING REMARKS 

In this paper, we have studied the effective spin- 
orbital-lattice coupled model which we derived to under- 
stand the peculiar phase transition to the spin-singlet 
state in ATiSi206 {A = Na, Li). Using the mean- field- 
type analysis and the numerical quantum transfer ma- 
trix method, we have clarified that the interplay between 
spin and orbital degrees of freedom plays a central role in 
thermodynamic properties of the system. At high tem- 
peratures, both spin and orbital correlations are anti- 
ferro type and compete with each other. As temperature 
decreases, the antiferromagnetic spin correlation grows 
rapidly and yields the sign change of the orbital corre- 
lation from antiferro to ferro type so that the frustra- 
tion is released. This ferro-type orbital correlation with 
the Jahn- Teller coupling finally causes a transition to the 
spin-dimer and orbital-ferro ordered phase. The feedback 
effect between orbital and spin degrees of freedom results 
in peculiar temperature dependence of the magnetic sus- 
ceptibility. We have shown that the magnetic susceptibil- 
ity data for NaTiSi206 can be explained by our effective 
model with realistic values of parameters. 

The transition to the spin-singlet state in the present 
system shows several different aspects from conventional 



spin-Peierls systems. One is the temperature dependence 
of the magnetic susceptibility. In the present system, the 
rapid decay of the susceptibility due to the spin-singlet 
formation occurs even without a broad peak as a fin- 
gerprint of well-developed spin correlations. This is be- 
cause the driving force of the transition is not the mag- 
netoelastic interaction but the orbital-ferro correlations 
assisted by the Jahn- Teller distortion. However, this does 
not mean that the orbital-lattice physics is dominant as 
in Cg electron systems such as CMR manganites. In 
the present <2g electron system, the Jahn- Teller energy 
scale is considerably smaller than the orbital and spin 
exchange interactions, and the orbital-ferro correlation 
is induced by the keen competition between the orbital 
and spin degrees of freedom. These illuminate a unique 
feature of the present t2g system, namely, there the inter- 
play between orbital and spin appears explicitly without 
being dominated by Jahn- Teller physics. 

Another peculiar aspect of the spin-dimer state in 
our model is that the spin-singlet pairs are formed on 
the longer Ti-Ti bonds rather than the shorter ones. 
In the present system, say, the dxy orbital ordering is 
concomitant with the flattening of TiOe octahedra in 
the z direction which elongates Ti-Ti bonds in the xy 
plane. Since the spin exchange interaction is strong be- 
tween the nearest-neighbor sites in the xy plane in this 
dxy ordered state, the spin-singlet dimers are formed on 
the longer Ti-Ti bonds. This aspect is clearly different 
from the conventional spin-Peierls systems in which the 
spin-singlet pairs are on shorter bonds. In our model, 
however, we take account of only the tetragonal Jahn- 
Teller mode which couples to dxy and dy^ orbitals. For 
detailed comparisons with the experimental data which 
show much complicated lattice structure at low temper- 
atures, it is necessary to include more general lattice dis- 
tortions in our theory. In particular, we note that the 
magnetoelastic coupling can be substantial in the low- 
temperature orbital-ordered phase since the orbital-ferro 
ordering largely enhances the effective spin exchange 
coupling on strong bonds up to ^ lOOOK as shown in 
Sec. IIV B 21 The magnetoelastic coupling will cause an 
opposite effect on the Ti-Ti bond lengths since it tends 
to shorten the spin-singlet Ti-Ti bonds. Further study is 
necessary to conclude the low-temperature lattice struc- 
ture. Note that the inclusion of such magnetoelastic ef- 
fect does not alter our conclusions on the mechanism of 
the phase transition because it becomes important only 
when the orbital-ferro ordering is well established far be- 
low the critical temperature. 

In the present study, we have compared our results of 
the magnetic susceptibility with the experimental data 
for the Na compound NaTiSi2 06. There is another com- 
pound in this pyroxene family, i.e., LiTiSi2 06. The Li 
compound also shows the phase transition at Tc = 230K 
showing a sudden decay of the magnetic susceptibility 
below Tc.^ However, the susceptibility data show some 
extra anomalies probably due to impurity phases. We 
believe that our model describes essential physics in both 
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the Na and Li compounds. Further experimental study 
including the sample refinement is desired to compare the 
data of Li compounds to our results. 

As shown in Sec. IIII Al and in Appendix B, the low- 
temperature phase of our effective model turns into the 
spin-ferro and orbital-antiferro ordered state for larger 
values of 77 than the critical value of rjc — 0.18. The pa- 
rameter 77 is the ratio of the Hund's-rule coupling to the 
intraorbital Coulomb repulsion, and is considered to be 
r] 0.1 < rjc in the present compounds. Unfortunately, 
it is difficult to control the parameter rj experimentally, 
however, it may vary for different compounds to some ex- 
tent. If rj becomes close to r]c in some compound, there 
would occur interesting phenomena related to the criti- 
cality of the phase transition at r]c. One interesting ex- 
ample is the phase transition from the spin-dimer and 
orbital-ferro state to the spin-ferro and orbital-antiferro 
state by applying the external magnetic field. We have 
investigated this issue in our effective model, and indeed 
observed the field-induced phase transition. The results 
will be reported elsewhere. The experimental study of 
this issue is left for future study. 
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APPENDIX A: QUANTUM TRANSFER MATRIX 
METHOD 



In this appendix, we briefly review the algorithm of 
the QTM method.3^ The QTM calculation is apphed to 
the effective ID model (|36|) in the self-consistent scheme 
as described in Sec. IIV Al 

After the Suzuki- Trotter decomposition, the partition 
function is represented in terms of the transfer matrix as 



L/2 



Z = Tt e"^^ = hm Tr TT Tm 



(Al) 



where /3 — 1/T (we set the Boltzmann constant — 1), 
M is the Trotter number, and L is the system size. The 
transfer matrix Tm is given by 



M 



Tm= e-'3''^"-i'^"/^^e-''''^-^"+i/^ , (A2) 

where Ti. is decomposed into the summation of the local 
Hamiltonian h^ i^i for bond (i, i + 1). We omitted the in- 
dex n on the transfer matrix Tm since T/\/ for the present 
system is invariant under the two-site translation. 



The advantage of the QTM method is that we can 
calculate thermodynamic quantities directly from the 
largest eigenvalue Amax and the corresponding right \v^) 
and left (w'l eigenvectors of T/\/.^^ (Note that the eigen- 
vectors !«'") and (w'l are different in general since Tm is 
an asymmetric matrix.) The partition function of the 

infinite system is found to be Z = lim^^joo Amax, and 
consequently, the free energy per site can be obtained as 



(A3) 



Hence, once the value of Amax is obtained, we can cal- 
culate any bulk quantity in the thermodynamic limit by 
taking the appropriate derivative of the free energy /. 

Furthermore, we can calculate expectation values of 
site and bond operatorsi^ For instance, an expectation 
value of a product of operators C'2n-i,2nC'2n,2n+i defined 
at sites (2n — 1, 2n, 2n -|- 1) is obtained from the formula 



(C2ri-l,2n02 n,2n+l / 



lim 

A/— >oo Amax 



(A4) 



where 



2n—l,2n 

o Tm- (A5) 



In the present study, we diagonalize the transfer matrix 
Tm numerically using the simple power method, which is 
known to be stable in the diagonalization of asymmetric 
matrixiSS Then, using Eqs. IIA4p and ljA5|l . we calculate 
the magnetization m [Eq. 139|) ] and the NN spin and 
orbital correlations Cg and Co [Eqs. and ||1T|) ] by 



(A6) 



^ _ {S2n-1 ■ S2n) + {S2n ■ S2n+l) 

— , 



{T2n-lT2n) + {T2nT2n+l) 



(A8) 



respectively. The NN spin correlations on odd and even 
bonds [Eqs. and (|^ ] are given by 



Cocid / c* \ 

s — \J2n-l ■ J2nl 

Ccvcn 1 c C \ 

s — \J2n ■ J2n+ll 



(A9) 
(AlO) 



respectively. In the spin-dimer phase, C°'^'^ and C° 
may take different values. 



APPENDIX B: TRANSITION TO SPIN-F AND 
ORBITAL-AF PHASE 

In this appendix, we discuss thermodynamic proper- 
ties of the model |^ for larger values of rj than those 
studied in Sec. IIV Bl There, the sF-oAF ground state is 
expected to be stable from the mean-field-type analysis 
in Sec lmH 
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FIG. 9: Temperature dependences of (a) the magnetic sus- 
ceptibility Xi (b) the nearest-neighbor spin correlation Cs, 
and (c) the nearest-neighbor orbital correlation Co for rj = 
0.16, 0.18 and 0.20. Symbols represent the results for M = 4 
while solid and dotted curves are those for M — 3 and M = 2, 
respectively. 



We first discuss our numerical results for the spin- 



orbital model T^so, i-e., the effective model (|36|) without 
the JT coupling 7 = 0, corresponding to the results in 
Sec. IIVB II Figure El (a) shows the results of the mag- 
netic susceptibility x- The result for 77 — 0.16 exhibits a 
sharp drop as T ^ similarly to the results in Fig.O in- 
dicating the spin-singlet ground state. On the contrary, 
for rj — 0.18 and 0.20, x exhibits a divergent behavior 
as T — > 0, suggesting that the ground state is mag- 
netic. These suggest that there is a ground-state phase 
transition between nonmagnetic and magnetic phases at 
r? = 7?c - 0.18. 

To clarify the nature of the transition in more detail, 
we show the results of NN spin and orbital correlations 
Cs and Co in Figs. El (b) and (c), respectively. The results 
for 77 — 0.16 show similar features to those for smaller rj 
shown in Sec. lIVBl and indicate that the system exhibits 
the sD-oF ground state. For rj = 0.20, on the other hand, 
the antiferro-type orbital correlation Co develops rapidly, 
yielding the sign change of the spin correlation Cg . After 
the sign change, the two correlations grow cooperatively, 
and approach the values of Cs = 1/4 and Co — —1/4 
as r — * 0. These are the values expected in the sF-oAF 
ground state shown in Fig. [21(b). Hence, the system with 
77 = 0.20 exhibits the sF-oAF ground state. For r] = 0.18, 
Cs and Co are both antiferro type and compete with each 
other down to the lowest temperature studied here. This 
suggests that 77 = 0.18 is close to the phase boundary 
between the sD-oF and sF-oAF phases. Therefore, we 
conclude that the ground state of the spin-orbital model 
Tiso undergoes a phase transition between the sD-oF and 
sF-oAF phases at rj — rjc ~ 0.18. The critical value 
is in good agreement with the mean-field prediction in 
Sec. im^ 

In the case of a finite JT coupling, we have found that 
there occurs a finite-temperature phase transition to the 
low-temperature sF-oAF phase for rj > rjc or to the sD-oF 
phase for rj < rjc- The details will be reported elsewhere. 
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